Jet Quenching and the p > n Anomaly at RHIC 
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PHENIX data on Au + Au at y/s = 130 AGeV suggest that p yields may exceed ir~ at high 
Pt > 2 GeV/c. We propose that jet quenching in central collisions suppresses the hard PQCD 
component of the spectra in central A + A reactions, thereby exposing a novel component of baryon 
dynamics that we attribute to (gluonic) baryon junctions. We predict that the observed p > n~ and 
the p > tt + anomaly at pt ~ 2 GeV/c is limited to a finite pt window that decreases with increasing 
impact parameter. 



PACS numbers: 12.38.Mh; 24.85.+P; 25.75.-q 

Introduction. Recent data on Au + Au reactions at 
y/s = 130 AGeV from the Relativistic Heavy Ion Col- 
lider (RHIC) have revealed a number of new phenomena 
at moderate high px ~ 2 — 6 GeV/c. The high px spectra 
of 7T° in central collisions were found by PHENIX [ to 
be suppressed by a factor ^3 — 4 relative to perturbative 
QCD (PQCD) predictions scaled by the Glauber binary 
collision density (Ta^d)). Also in non-central collisions, 
STAR [ | found that the linearly increasing azimuthal 
asymmetry moment, U2(px)> saturates at ~ 0.15 forpx > 
2 GeV/c in contrast to predictions based on ideal hydro- 
dynamics [ ||. From Rcf. [ ^] these features were in- 
terpreted as evidence for jet quenching in a dense gluon 
plasma with rapidity density dN 9 jdy ~ 1000. The initial 
proper gluon density may thus have reached far into the 
deconfined QCD phase with p g (r ~ 0.2 fm/c) > 10/fm 3 . 
Jet quenching at RHIC, if confirmed by further measure- 
ments, opens the door to a new class of diagnostic tools 
that probe AA dynamics and the transient quark-gluon 
plasma created in such collisions. 

One of the unexpected results reported by PHENIX 
and STAR is that in contrast to the strong ir° quench- 
ing at pt > 2 GeV/c, the summed negatively charged 
(tt~ + K~ + p) and the corresponding positively charged 
hadron spectra were found to be quenched by only a fac- 
tor ~ 2 [ [I], 9. Even more surprisingly the identified 
particle spectra analysis at PHENIX [ g| suggests that 
Rb(pt) = _p/7T~ > 1 at pt > 2 GeV/c. Thus, baryon and 
anti-baryon production may in fact dominate the mod- 
erate high pt hadron flavor yields, a phenomenon never 
before observed. 

These and other data point to novel baryon transport 
dynamics playing role in nucleus- nucleus (A A) reactions. 
An important indicator of this are STAR [ [jj data that re- 
vealed a high valence proton rapidity density (~ 10), five 
units from the fragmentation regions, and a p/p ~ 0.65 
at mid-rapidity. An attractive dynamical model that ex- 
plains copious mid-rapidity baryon and anti-baryon pro- 
duction is based on the existence of topological gluon field 
configurations (baryon-junctions) [ |j, ^|, |l(|. Junctions 
predict long range baryon number transport in rapidity 
as well as hyperon enhancement (including f2 _ ) [ M and 
considerable pt enhancement [ |^, |l0| relative to conven- 
tional diquark-quark string fragmentation [ [ll]] . In this 



paper we propose that the baryon/meson anomaly is due 
to the interplay between the jet quenched [ [|, O, |lj] 
hard component (dNh) and the phenomenological soft to 
moderate px component (dN s ) of p and tt~ . The ob- 
served baryon junction component is assumed to exist in 
nucleon-nucleon (NN) and well as AA reactions, but in 
NN [ ||, |^| it has a smaller amplitude, and its contribu- 
tion to the high px baryon spectrum is "obscured" by 
unquenched mini-jet fragmentation into pions. 

In this paper we extend the study of K^/tt^ 1 [ 
and p/p [ |l4| ratios to compute the expected differential 
baryon flavor ratio, Rb(pt) = P/^~ ■ We also generalize 
the two component soft+hard dynamical model of [ Q to 
include the novel baryon junction component. We pre- 
dict that the enhancement of the p/tt~ ratio reported by 
PHENIX is actually limited to a finite moderate px range 
2-6 GeV/c. Beyond this range, the hadron ratios are ex- 
pected to converge back to the quenched PQCD base. In 
Au + Au, our proposed mechanism can be further tested 
via a predicted systematic decrease of p/tt~ with increas- 
ing impact parameter, i.e. decreasing nucleon participant 
number. 

Reference PQCD flavor composition. The stan- 
dard PQCD approach expresses the differential hadron 
cross section in NN — > hX as a convolution of the mea- 
sured structure functions f a /N{xa,Q 2 a ) for the interact- 
ing partons (a = a, b), with the fragmentation func- 
tion D h / c (z,Qc) for the leading scattered parton c into 
a hadron of flavor h and the elementary parton-parton 
cross sections do-^ ah ^ cd ^ jdi. Folding over an intrinsic par- 
tons (Gaussian) /(kx) distribution: 
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-S(s + u + t), (1) 



where x a , x b are the initial momentum fractions carried 
by the interacting partons, z c — Ph/p c is the momen- 
tum fraction carried by the observed hadron. For NN 
we use phenomenological smearing (k T ) = 1.8 GeV 2 /c 2 . 
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Comparison between the PQCD calculation ([y) and the 
charged hadron multiplicities measured by the ISR (pp) 
and UA1 (pp) experiments [ [l5| at y/s = 53 GeV and 
200, 900 GeV respectively are in good agreement [ [To| for 
various structure and fragmentation function parameter- 
izations [ [u| 0. 
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FIG. 1. The ratio Rb ~ (p+p)/ ( 7r+ + t ) is shown for pp 
at y/s = 1.8 TeV as a function of p T . Low p T FNAL E735 
data [ [if] is shown for comparison. String fragmentation of 
mini-jets via HIJING1.35 [ |ll| is also shown. 

Fig. 1 compares Rb(pt) from the conventional inde- 
pendent and string fragmentation PQCD phenomenolo- 
gies to available E735 [ [if] data from yfs = 1.8 TeV 
pp collisions. While the pion dominated charged particle 
inclusive data are fit well by ([!]) with different fragmenta- 
tion models [ 
predicted p/ir 
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there is a large difference between the 
ratios in the shown px range. Comparing 
results from two different parametrizations of fragmenta- 
tion functions [ [l7| and the HIJING1.35 string fragmen- 



tation event generator [ 11 reveals a rather large theo- 
retical uncertainty ~ 3 on Rb at high px. Nevertheless, 
the prediction that Rb $ 0.5 at high px is robust and in- 
tuitively reasonable given that both quark and gluon jets 
prefer to fragment into lighter mesons than baryons. We 
conclude that a systematic study of the px dependence of 
the baryon/meson ratios at RHIC should include the ~ 3 
uncertainty of the baryon fragmentation schemes shown 
by an error band in Fig. 1, but the PQCD predictions for 
P/tt~ is that it stays below 0.5 for px > 2 GeV/c. 

Anomalous baryon component. We generalize 
next the (soft+hard) two component model in Ref. [ Q| 
and test it against the observed charged particle quench- 
ing pattern [|,|§. 

In AA collisions the semi-hard PQCD component is 
reduced due to the non-abelian energy loss of hard par- 
tons before fragmentation as in [ [| . In the calculations 
below, we use radiative spectrum dl/dx (corresponding 
to fractional energy loss AE/E) computed in the GLV 
(Gyulassy-Levai-Vitev) formalism [ The numerical 
techniques developed in [ n2| allow us to extend the jet 



quenching computation to partonic energies as low as 
~5 GeV that are relevant for the px range addressed in 
this paper. Jet production and propagation through the 
medium is treated as in [ Q including the realistic Wood- 
Saxon nuclear geometry. We also include 1+1D Bjorken 
and 1+3D Bjorken+transverse expansion of the interac- 
tion region. In [ [l^] the azimuthally averaged energy loss 
was found both analytically and numerically to have little 
sensitivity to transverse flow 

(AE 3D )^ (A£ 1D ) . (2) 

In the computation of the PQCD component we take 
into account nuclear shadowing as in [ [ll| and Cronin 
effect as in Ref. [ ^0|. In addition, we take into account 
multi-gluon fluctuations of the energy loss via the opacity 
renormalization approach as discussed in [ . Fluctua- 
tions in the GLV approach effectively reduce the opacity 
X oc dN 9 /dy by a factor ~ 0.5. The fitted initial gluon 
density quoted in Fig. 2 below takes that reduced energy 
loss effect into account. 

We parameterize the soft phenomenological component 
of moderate px hadrons as follows: 



dN s (b) __ ^ dn a e -PT/T=(b) 



dyd 2 pq 



(3) 



As in string models the soft component is assumed to 
scale with the number of participants (N part ). In Eq. (|f) 
we also account for the possibly different mean inverse 
slopes T a for baryons and mesons. In the junction pic- 
ture [ H , the large T p may arise from the predicted [ [|] 
smaller junction trajectory slope a'j ps ct' R /3. This im- 
plies that the effective string tension is three times higher 
than l/(2Tra' R ) ~ 1 GeV/fm leading in the massless limit 
to (p\)j ~ 3 (Px)i?. In terms of the string model the 
factor three enhancement of the mean square px is due 
to the random walk in px arising from the decay of the 
three strings attached to the junction. Naively, we would 
thus expect T p ~ V3T 71 ". In a detailed Monte Carlo 
study [ using HIJING/BB [ |9| this relation was in- 
deed found to hold in the 1 < px < 2 GeV/c range with 
T" ps 220 MeV and T p ps 370 MeV - VST* with junc- 
tions included and T p ~ 270 MeV without. At higher 
px the PQCD mini-jets cause the apparent pion slope to 
increase systematically with px, while at lower px the 
resonances cause T n to be smaller. In going to periph- 
eral collisions a small decrease in the mean inverse slopes 
AT 7r- = _ 1Q MeV and AT p _ _ 5Q MeV for approxi _ 

mate consistency with pp data [ ^] is introduced. 

The general hydrodynamic solution [ ^3) including 
transverse flow is of course much more complex at low 
px than our simplified parameterization (^). However, at 
large px all hadronic inverse transverse slopes tend to a 
Doppler shifted value (T e ff = Tt exp(r/ r ), in terms of the 
transverse flow rapidity r\ r at a freeze-out isotherm Tf). 
As we show below, boosted thermal sources including rel- 
ativistic transverse flow does not predict an anti-baryon 
anomaly in the px ~ 2 GeV/c range [ 10 . 
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The baryon transport mechanism suggested in [ P pre- 
dicts that the valance baryon number per unit rapidity 
has the form 

dN B dNP dNP nrr cosh(l-a B (0))y 

= P Z 77^? ) ( 4 ) 



dy 



dy dy 



sinh(l-a B (0))r m 



where a_e(0) — 1/2 is the baryon junction exchange inter- 
cept. For Au + Au reactions at y/s = 130 AGeV Eq. (@) 
predicts dN B Jdy ~ 10 which is in good agreement with 
RHIC data [fjj. 

Our phenomenological soft string and baryon junction 
components depends on 2 parameters - the total charged 
particle rapidity density and the mean inverse slope for 
pions. We take dN ch /dy ~ 650 for the top 6% central 
events from experiment [ [2f| and T" ~ 220 MeV. It 
has been argued that the soft string and baryon junction 
dynamics is manifested in nuclcon-nuclcon collisions [ 

li- 
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Therefore the existing pp data 
us with reference integrated particle ratios and predic- 
tions (e.g. p/ir~ — 0.06) for the case of AA reactions. 
These predictions are found to be consistent with the 
measurements at RHIC within 30%. Deviations of those 
Pt integrated ratios from the pp case are expected (e.g. 
strangeness enhancement) but they cannot account for 
the factor of ~ 5 anomalous enhancement of the pt dif- 
ferential baryon/meson ratio reported by PHENIX [ [|. 

The full hadronic spectrum is a sum of soft and hard 
components 



dN{b) _ dN s (b) dN h (b) 
dyd 2 pT dyd 2 pt dyd 2 pT 



(5) 



with a hard component computed from the cross section 
in nucleon-nucleon collisions modified by the medium in- 
duced energy loss [ |12, 13 and the Glauber profile den- 



sity at a given impact parameter. To avoid over-counting, 
we smoothly turn on the power law PQCD component at 
p T ~ 2 GeV/c [|. 

We first compute the ratio Raa(pt) of the hadron mul- 
tiplicity in Au + Au at y/s = 130 AGeV normalized to 
the corresponding multiplicity in pp reactions scaled by 
the number of binary collisions. Such scaling is relevant 
for the comparison of the moderate to high pt parts of 
the hadronic spectra and in absence of nuclear effects 
is expected to be unity. Fig. 2 shows Raa m the 10% 
central collisions for neutral pions and inclusive charged 
hadrons. We note that including only nuclear shadow- 
ing and Cronin effect gives a ratio consistent with unity 
within 10% for pt > 2 GeV/c and inconsistent with data 
if jet quenching is not taken into account [ ^?J. The fit- 
ted parton energy loss is specified by initial gluon rapidity 
density dN 9 /dy = 800 ± 100 and this density is close to 
predictions based on saturation models [ |2^| . This high 
gluon density takes into account the reduction of the ef- 
fect of induced radiative energy loss due to multi-gluon 
fluctuations as discussed in detail in Ref. [ ^l| . 

In relativistic heavy ion collisions early thcrmalization 
of the soft background partons is expected. This allows 



the estimate of the Debye screening /j, (which is a natu- 
ral infrared cut-off in medium) from thermal pQCD. The 
screening scale then enters the partonic cross section and 
the transport coefficient used in energy loss calculations [ 
|l2| . Larger variations of the initial density of the medium 
have also been studied [ |l(| but they give a less satisfac- 
tory description of the PHENIX data [ |IJ . 

At pt > 2 — 3 GeV/c pions are perturbative and ex- 
hibit an almost constant suppression as predicted by the 
GLV formalism [ [l2| [ll| . In contrast at moderate high pt 
the non-perturbative baryon junction component is unaf- 
fected. Thus in central collisions baryon excess at inter- 
mediate transverse momenta is expected relative to NN. 
This results in a smaller degree of suppression of inclusive 
charged hadrons which is pt dependent. At larger trans- 
verse momenta baryon production is expected to be dom- 
inated again by PQCD (as in Fig.l). Therefore baryons 
are expected to show comparable suppression to pions at 
high pt > 5 GeV/c. From Fig. 2 we conclude that such 
physical picture is compatible with the data and can be 
easily tested with future data in the pt > 5 GeV/c range. 
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FIG. 2. The ratio of charged hadron and it multiplicities 
to the binary collision scaled pp result is shown from [ 
The curves utilize the GLV quenched hard spectrum and 
the modified soft component Eq. (^). 

The p/tt~ anomaly. In our discussion of the pt de- 
pendence of p/V - ratio at RHIC we first turn to a simple 
description based on boosted thermal sources 

dN ~ mvKi {^^t) lo , (6) 

where tanh?7 r = v±. A computation of Rb{pt) with 
Tf = 160 MeV and v± = 0.6 is included in Fig. 3. It pre- 
dicts p/n~ < 1 in the intermediate pt window and grows 
monotonically with pt to Rb = 2 which is a general fea- 
ture of hydrodynamic calculations. We note however, 
that variations of hydrodynamic initial density profiles, 
equation of state, and freeze-out criteria as well as dissi- 
pative effects may also play a role in the p/tt~ anomaly, 
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since Rb < 1 may also occur in hydrodynamics with 
late chemical freeze-out [ || as well as in the hybrid hy- 
dro+UrQMD model due to dissipative effects [ |24| . 

As a further test of our model we consider next the 
expected centrality dependence of Rb- The p/-K~ ra- 
tio computed for 3 different centralities in the soft+hard 
model Eq. (|^) is also shown in Fig. 3. In central collisions 
the interplay between the anomalous baryon component 
of p and the quenched PQCD component of ir~ leads 
to maximum of Rb near pt ~ 3 — 4 GeV/c. At large 
Pt > 5 — 6 GeV/c we predict a gradual decrease oi Rb 
below unity consistent with the the PQCD baseline cal- 
culations (see Fig. 1). In Fig. 3 we have included through 
error bands the factor of ~ 3 uncertainty in the fragmen- 
tation functions into p at high pt. The solid and dashed 
curves reflect the difference between the N part and 
scaling of the junction component 

In peripheral reactions the size of the interaction re- 
gion as well as the initial density of the medium de- 
crease leading to reduction of energy loss. The absence 
of quenching reduces the observability of the anomalous 
component and the p/7r~ ratio may stay below unity for 
all pt- The case of peripheral reactions is hence similar 
to pp collisions. The experimentally testable prediction 
of the model (^|) is therefore that the maximum of the 
Rb = p/vr - ratio decreases with increasing impact pa- 
rameter, decreasing participant number, or equivalently 
decreasing dN ch /dy. 
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FIG. 3. The centrality dependence of Rb(pt) is predicted 
for three different centralities. Solid (dashed) lines corre- 
spond to A 1 (A A / 3 ) scaling of the junction component. The 
ratio of p and ir~ fits to PHENIX data on central reactions 
is shown for comparison. A boosted thermal source (dashed 
line) is also shown. 

Conclusions. We found that the interplay between 
jet quenching at RHIC computed as in [ Q using the 
GVW, GLV formalisms [ [|, |l2| and a postulated novel 
baryon junction component could account for the pt > 
2 GeV/c p/ir~ £ 1 puzzle suggested by the data. We pro- 
pose that in central collisions jet quenching exposes this 



new baryon junction component [ g, ||. This was shown 
in Fig. 2 to provide a natural explanation for the smaller 
effective quenching of moderate pt charged particles than 
7T°. Fig. 3 shows that Rb could reach values well above 
unity in a finite moderate pt domain but is predicted to 
decrease to well below unity at pt > 5 GeV/c in contrast 
to thermal or hydrodynamic models. Whereas the maxi- 
mum baryon excess, Rb max, depends on unknown details 
of the baryon junction component the shape of Rb(pt) 
and its specified centrality and transverse momentum de- 
pendence are largely a consequence of the jet quenching 
phenomenon similar to the saturation of V2(pt) reported 
by STAR [§. 

Since from pt integrated hadron yields p/p ~ 0.65 [ 
[l], |^], the predicted pt dependence of baryon enhance- 
ment and subsequent reduction at high pt may in fact 
be more readily observable in the p/n + ratio. 
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FIG. 4. The centrality and pt dependence of the p/n + ratio 
is predicted for the same centrality classes as in Fig. 3. The 
solid (dashed) lines compare the difference between N part 
(Npirt) scaling of the junction component and 3(1) x the 
quenched fragmentation component as in Fig. 3. 

We find that the p/ir + ratio is typically ~ 50% larger 
than p/ir~ as a function of pt at moderate transverse 
momenta. Comparison between Fig. 3 and Fig. 4 shows 
that for pt > 6 GeV/c the baryon/meson ratio for posi- 
tive hadrons predicted here is ~ 2 — 3 times bigger than 
the corresponding ratio for negative hadrons. This is due 
to the valance quark dominance in the PQCD fragmen- 
tation into baryons versus antibaryons at high px (see 
also Refs. [ [l0| [l4|]). In contrast, hydrodynamic mod- 
els with flavor independent freeze-out would predict an 
asymptotic ratio ~ 2. 

In summary, a novel test of the baryon junction hy- 
pothesis [ |^, |i| for baryon number production and trans- 
port in ultra-relativistic nuclear collisions has been pro- 
posed through the predicted anomalous dependence of 
high px p/tt - and p/ir + ratios and their centrality de- 
pendence. The existing PHENIX data [ ||] strongly mo- 
tivate the extra effort needed to measure (anti)baryon 
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spectra into the pr > 5 GeV/c range as a complement 
to the (anti)baryon number rapidity transport measure- 
ments by STAR [ § and soon BRAHMS [ §§. 
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